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Abstract 

We analyze a sample of 3 million quantum-correlated D°D° pairs from 818 pb _1 of e + e~ collision 
data collected with the CLEO-c detector at E cra = 3.77 GeV, to give an updated measurement 
of cos 5 and a first determination of sin 5, where S is the relative strong phase between doubly 
Cabibbo-suppressed D° — > K + ir~ and Cabibbo-favored D° — > K + ir~ decay amplitudes. With no 
inputs from other experiments, we find cos 5 = 0.81^1^^1, sin 5 = —0.01 ± 0.41 ± 0.04, and 
|<5| = (10^53^J 3 )°. By including external measurements of mixing parameters, we find alternative 
values of cos 5 = 1.15±g;i?±g:gg, sin<5 = 0.56^;^;^, and 5 = (18±^)°. Our results can be used 
to improve the world average uncertainty on the mixing parameter y by approximately 10%. 
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t Now at: National Taiwan University, Taipei, Taiwan 
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I. INTRODUCTION 



Charm mixing in the Standard Model is conventionally described by two small dimen- 
sionless parameters: 

x = 2 — ^ (1) 

r 2 — Ti 

y = r t r . ( 2 ) 
1 2 + 1 1 

where M 12 and r 12 are the masses and widths, respectively, of the neutral P meson CP 
eigenstates, Pi (CP-odd) and P 2 (CP-even), defined by 

IP ) + IP ) 

W - (3) 

IP ) - IP ) 

W = ^LJ, (4) 

assuming CP conservation. The mixing probability is then denoted by Pm = (x 2 + y 2 )/2, 
and the width of the P° and P° flavor eigenstates is T = (I\ + r 2 )/2. 

Recent experimental studies of charm mixing parameters have probed x and y directly [T]- 
[7], as well as the "rotated" parameter y' = ycos5 — x sin S [SHE]. Here, —5 is the relative 
phase between the doubly Cabibbo-suppressed P° — > K + tt~ amplitude and the correspond- 
ing Cabibbo-favored P° -> K + ic~ amplitude: (A + 7r _ |P°)/(A" + 7r~|P°) = re~ iS . We adopt a 
convention in which 5 corresponds to a strong phase, which vanishes in the SU(3) limit [TT] . 
The magnitude r of the amplitude ratio is approximately 0.06. In this article, we update an 
analysis [12] that first directly determined cos 5 using correlated production of P° and P° 
mesons in e + e~ collisions produced at the Cornell Electron Storage Ring and collected with 
the CLEO-c detector. In the current analysis, we also present a first measurement of sin S. 

At the ^(3770) resonance, the P°P° pair is produced with no accompanying particles, 
so it is in a quantum-coherent C = — 1 state. As a result, the P°P° decay rates differ from 
incoherent decay rates because of interference effects. These differences depend on y (to 
first order) and on strong phases of the decay amplitudes [HI IT3H22] . As in our previous 
analysis, we implement the double-tagging method with a \ 2 fit described in Ref. |23j, 
where, in addition to extracting the number of P°P° pairs produced (AT) and the branching 
fractions (£>) of the reconstructed P° final states, we simultaneously determine y, x 2 , r 2 , 
cos S, and sin S, all without needing to know the integrated luminosity or P°P° production 
cross section. The main improvements in the current analysis are: use of the full CLEO-c 
^(3770) dataset (which is three times larger than the previous dataset), reconstruction of 
additional CP eigenstates, reconstruction of semimuonic P° decays, addition of modes that 
provide sensitivity to sin 5, and direct measurement of the amplitude ratio r. 

As before, we neglect CP violation in P decays and mixing. Recently, evidence has 
been found [21H2S] for direct CP violation in P — > K + K~ and P — > 7r + 7r~ decays, with 
CP asymmetries of O(10~ 2 ). The current analysis also uses P — > K + K~ and P — > 7r + 7r~ 
decays, and the above CP asymmetries would bias our determinations of cos 5 and sin 5 
from these modes by 0(1O~ 2 ). By itself, this bias is much smaller than our experimental 
uncertainties, but its effect is further diluted by our use of additional CP eigenstates in the 
analysis. CP violation in mixing and in the interference between mixing and decay would 
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bias our measured value of y. However, we are insensitive to these sources of CP violation 
at the levels currently allowed by experimental constraints [27]. 

The paper is organized as f ollow s. In Section [TT| we review the formalism of quantum- 
correlated D°D° decay. Section III describes the event selection criteria and D reconstruction 
procedures. The external measurements used in the fit are summarized in Section IV 



Systematic uncertainties, which are also input to the fit, are discussed in Section N\ Finally, 



we present and discuss our main fit results in Section VI 



II. FORMALISM AND PHENOMENOLOGY 

For decays of isolated D° mesons, we define the following quantities for each final state i: 



r 2 = 



JA(x)^(x)rfx 
/^(x)^(x)dx 
/ ,l,(x),l;(x)r/x 



(5) 



1 -rJA 4 (x)A*(x)dx W 

where A{ = (i\D°) and A% = (i\D°) are the amplitudes for the transitions of D° and D°, 
respectively, to the final state i. The integrals are taken over the phase space for mode i. 
Thus, 5i is an average phase for the final state i, and Ri G [0, 1] is a coherence factor [28] 
that characterizes the variation of <5« over phase space. If the final state is two-body, like 
K~tt + , then 5, is constant over phase space, and R4 — I. 

For a D°D° pair produced through the ^(3770) resonance, the decay rate to an exclusive 
final state {i, j}, where i and j label the final states of the two D mesons, follows from the 
antisymmetric amplitude M.%f 

T{i,j)<xM% = lAiAj-AiAjf 

= KWOW " <«|A>0W| a + 0(x 2 ,y 2 ), (7) 

where the 0(x 2 , y 2 ) term represents a mixed amplitude. The interference between mixed and 
unmixed amplitudes vanishes when integrated over time because it depends on the difference 
between the D° and D° decay times. If we denote the charge conjugates of modes % and j 
by 1 and J, then Equation ^ leads to the following expressions in terms of the parameters 
defined above: 

T(i,j) = T(t,j) oc A 2 A 2 (l + r 2 r 2 — 2rji?j cos 5i rjRj cos 5j — IriRi sin 8i rjRj sin 5jj (8) 
T(i,j) = oc A 2 A 2 (r 2 +r| — IriRi cos5j rjRj cos5j + 2r{Ri sin5j rjRj sin 5 A , (9) 

where the latter rate is reduced by half if i and j are identical. Experimentally, these rates 
correspond to yields of double tags (DT), which are events where both D° and D° are 
reconstructed. 

The above amplitudes are normalized to the uncorrelated branching fractions Bf. 

B, = B(D° i) = A 2 [1 + nRiiy cos Si + x sin fy] (10) 
Bf = B(D° i) = A 2 \r 2 + nRiiy cos 6i -xamSijl . (11) 
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These Bi are related to rates of single tags (ST), or individually reconstructed D° or D° 
candidates, which are obtained by summing over DT rates: 



T(i,X) = ]T + r(«,5)] =B i + B i = A 2 (l + 2yr i B 4C os8 i + r?) 



(12) 



Here, we have used an expression for y in terms of r», Ri, and 5i, which is derived from 
Eqs. pP) and Eqs. dlpb: 





'/ IA(x) 


- ^(x)| 2 dx - / |A(x) + A(x)| 2 dx 




JiA(x) 


- ^(x)| 2 (ix + / |A(x) + v4;(x)| 2 (ix 



' £^ 2 (i + rf) ' 



(13) 



Thus, both ST rates and the total rate, T^o^o, are unaffected by quantum correlations 
between the D° and D° decays, and our sensitivity to mixing comes from comparing ST to 
DT rates. 

Table H] gives the notation for the various r« and <5, that appear in this analysis. The final 
states of mixed CP (denoted by / and / below) that we consider are K^ir^ and K^ir+n^ . 
Following Ref. [29], the Kg7T + ii^ Dalitz plot is divided into eight bins according to the strong 
phase of the decay amplitude. We denote the portions of KgiT + Tr~ in phase bin i by Yi and 
Y~i, where m K o^- < m K o slT + for Y iy and m K o n - > m K o s ^+ for %. The corresponding amplitude 
ratio magnitudes and branching fraction ratios integrated over bin i are denoted by pi and 
Qi, respectively. As in Ref. J2H], we denote the real and imaginary parts of Eq. ^ by q 
and Si, respectively, but with the opposite sign convention for Sj. Semileptonic final states 
(£ ± ), CP-even eigenstates (S+), and CP-odd eigenstates (S-) have known values of and 
5i, which give them unique leverage in determining the parameters in the other final states, 



as demonstrated below. Note that, as shown in Eqs. 
equal rf in general. 



T0|) and pT]), the ratio Bj/Bi does not 



TABLE I. Parameters describing the ratio of amplitudes Ai and Ai for the final states i. The • 
indicates that we do not make explicit reference to 5i for the Y^ modes in this article, but consider 
only Ck and Sk instead. 



Final State Ti 5{ Ri cos Si Ri sin 5i BjjBi 





r 


5 


cos 5 


sin 5 




Y k /% 


Pk 




Ck 


Sk 


Qk 


S + 


1 


7T 


-1 





l 


S- 


1 





+1 





l 


e± 















Using the definitions in Table |l| we evaluate in Table [II] the quantum-correlated D°D° 
branching fractions, J-" cor , for all categories of final states reconstructed in this analysis; we 
also give the corresponding uncorrelated branching fractions, J-" unc . Comparing J 7001 with 
jrunc a n; ows us to extract y, r 2 , cos 5, and sin 5. Although we neglect x 2 and y 2 terms in 
general, we report a result for x 2 as determined solely from the suppressed {K ± n T , i^ ± 7r T } 
final states. 

From Table [TlJ one finds that, given r 2 and y, cos 5 can be determined by measuring 
the size of the interference between K~n + and a CP eigenstate. The CP of the eigenstate 
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TABLE II. Correlated (C-odd) and uncorrelated effective D°D° branching fractions, J- cor and 
J 711110 , to leading order in x, y and i?ws> divided by B{ for ST modes i and BiBj for DT modes 
Charge conjugate modes are implied. 

Mode Correlated Uncorrelated 

K--K+ 1 + i?ws 1 + #ws 

S+ 2 2 

S- 2 2 

n 1 + Q fc 1 + Q fe 

if - 7T+, i-C~7r+ i?M[(l + Rws) 2 - 4r cos(5(r cos 5 + y)] i? W s 

K _ 7r + , K + ir~ (1 + ^ws) 2 - 4rcos<5(rcos5 + y) 1 + ^ws 

K"vr+, 5+ 1 + i?ws + 2r cos 5 + y 1 + R ws 

R-TT+, S- 1 + i?ws - 2r cos <5 - y 1 + i? W s 

K~ir + ,£~ 1 — ry cos S — rx sin 5 1 

K~ir + , l + r 2 (1 — ry cos 5 — rx sin S) it/yvs 

K-7T+, (l+i?ws)(l + ft)-r 2 -^ l + i? ws Q. 
-2(r cos 5 + y)(piCi + y) + 2r sin 5 piSi 

K-ir+,Yi (l + Rws)(l + Q t )-l-r 2 p 2 R ws + Q, 
-2(r cos 5 + y)(piCi + y) -2r sin 5 piSi 

S+, S+ 1 

S-, S- 1 

S+, S- 4 2 

5+,r l + y 1 

S-,tr l-y 1 

5+, r» l + Q i + 2p i c i + y l + Qi 

S-,Yi l + Qi- 2pia - y l + Qi 

Yi, f 1 - piyci - piXSi 1 

Yi, i + p 2 (l - piyci - piXSi) Qi 



id) l J i Vj 

-2(piCi + y)(pjCj + y) + 2p i s i p j s j 
(l + Q^il + Q^-l-pfp 2 

-2{ Pi Ci + y)(pjCj + y)- 2p i s i p j s j 



V. V. ^T W T^j; ^ HiHj /0._LO- 



tags the K~tt + parent D to be a CP eigenstate with the opposite eigenvalue. Since this 
D eigenstate is a linear combination of the flavor eigenstates D° and D°, the decay rate is 
modulated by the relative phase between the D° — > K~n + and D° -+ K~ir + amplitudes. 

Similarly, probing sin<5 requires the interference of K~tt + with another mode, such as 
KgTT + ir~ , that has non-zero Ri sin Si. However, unlike CP eigenstates, the phases in Kg-K + 7T~ 
are not fixed by a fundamental symmetry, so we must measure sin S and Sj simultaneously. 
Since these sine factors only appear in products with other sine factors, there is an overall 
sign ambiguity, which can be resolved by combining our measurements of sin S and cos S 
with external measurements of y' and x. 

Our main source of information on y comes from CP-tagged semileptonic decays. In these 
weak transitions, the semileptonic decay width is independent of the parent D meson's CP 
eigenvalue. In contrast, the total width of the parent meson reflects its CP eigenvalue: 
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Ti,2 = r(l =F v)i so the semileptonic branching fraction for D\ or D2 is modified by a factor 
of 1 ± y. Thus, we determine y using exclusive final states {.S^^}, where the S± identifies 
the CP eigenvalue of the semileptonic decay's parent D. In this case, summing £ + and l~ 
rates gives F^e, ~ 2i3s ± i^(l ± y). 

For r 2 , we use the fact that, because of the vanishing interference between mixed and 
unmixed amplitudes, a DT with a semileptonic Kivn decay probes the bare matrix element 
squared, not the branching fraction, of the partner D. Therefore, we determine r 2 directly 
from {Kit, K£h>i} DT modes by taking the yield ratio for combinations with same-sign kaons 
and opposite-sign kaons. 



III. EVENT SELECTION AND RECONSTRUCTION 



Our current analysis is based on the full CLEO-c ^(3770) dataset with an integrated lumi- 
nosity of 818 pb -1 , collected with the CLEO-c detector, which is described in Refs. [30ti3l] . 
We estimate signal efficiencies, background contributions, and probabilities for misrecon- 
structing a produced signal decay in a different signal mode (crossfeed) using a GEANT- 
based (35] Monte Carlo simulated sample of uncorrelated D°D° decays with an effective 
integrated luminosity 20 times larger than that of our data sample. We reconstruct the 



7T + 7T 7T°. 



final states shown in Table III, with it — > 77, r] — > 77, K$ — > 7r + 7r~, and oj — > 
Final states without mesons and neutrinos are fully reconstructed. For modes with K\ 
mesons and neutrinos, which generally do not interact with the detector, we use a partial 
reconstruction technique, inferring the presence of the undetected particle via conservation 
of energy and momentum. In specifying the CP eigenvalue of a final state, we neglect CP 
violation in K° decays. 



TABLE III. D final states reconstructed in this analysis. 

Type Reconstruction Final States 

~T full K-7T+, Y -Y 7 

f full K+ir-,Y -Y 7 

S+ full K + K~, TT+TT-, K° S TT°TT 

S + partial A£tt°, #£77, K° L u 

S- full K° s tt°, K° sV , K° s lo 

S- partial A£vr 7r° 

£ + partial K~e + v e , K~ [i^v^ 

t~ partial A + e~P e , K+^v^ 



Final states that are common to those used in Ref. [12] are reconstructed with the same 
methods and selection criteria, except where noted below. In particular, the selection of 
7r , K , and K® candidates remains unchanged. For 7r° and rj candidates, we loosen the 
shower shape requirements to improve the agreement between efficiencies in data and those 
in simulated events. For all modes with oj candidates, we now apply a sideband subtrac- 
tion in the M(7r + 7r~7r°) spectrum. Figure [l] shows the invariant mass distribution of uj 
candidates, along with the signal region of 760.0 MeV/c 2 < M(7r + 7r~7r°) < 805.0 MeV/c 2 
and sideband regions of 600.0 MeV/c 2 < M(7r + 7r-7r°) < 730.0 MeV/c 2 and 830.0 MeV/c 2 < 
M(7r + 7r~7T°) < 852.5 MeV/c 2 . The limited range of the upper sideband is chosen to minimize 
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the effect of p° — > -k + -k~ and p — > 7r ± 7r° decays, which alter the shape of the background for 
M(7r + 7T~7r°) greater than approximately 870 MeV/c 2 . The sidebands are scaled by a factor 
determined by fitting the M(7r + 7T~7r°) distribution in simulated events to a signal Gaussian 
plus a polynomial background and integrating the fitted background function. This sideband 
subtraction eliminates peaking backgrounds, which accounted for 5 — 10% of the observed 
oo yields in Ref. |12] . We also make use of K^-k + -k^ DT yields, efficiencies, and background 
estimates from Ref. [29], for the subset of modes in that analysis without K l ti + tt~ , K~7t + tt 
or K~n + -K~it + . 
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FIG. 1. Distribution of M(7r + 7r~7r°) for to candidates. Data are shown as points with error bars, 
and the dashed lines mark the boundaries of the signal and sideband regions. As indicated in the 
text, the lower sideband ends at the lower limit of the graph. 



A. Single Tags 



We reconstruct ST candidates for the 8 modes in Table [TV] utilizing the technique de- 
scribed in Ref. [12j. We do not include ST yields for Kg7i + 7i~. As before, we identify ST 
candidates using two kinematic variables: the beam-constrained candidate mass M and the 
energy difference AE, which are defined to be 



M 
AE 



EZ/c* - 
Ed — E Q , 



Pl/c 2 



(14) 
(15) 
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where pd and Ed are the total momentum and energy of the D candidate, and Eq is the 



beam energy. After applying the mode-dependent requirements on AE listed in Table IV 



we determine the ST yields by fitting the M distributions, shown in Fig. [2j to a signal shape 
derived from simulated signal events and to a background ARGUS function [36J. 



TABLE IV. Requirements on AE for ST D candidates. 



Mode 


Requirement (GeV) 


K--K+ 


\AE\ < 0.0294 


K+TT- 


\AE\ < 0.0294 


K+K- 


\AE\ < 0.0200 


7T + 7r~ 


\AE\ < 0.0300 


K° s ir°ir 


-0.0550 < AE < 0.0450 


K° s tt° 


-0.0710 < AE < 0.0450 


K° S V 


-0.0550 < AE < 0.0350 




\AE\ < 0.0250 



The measured ST yields and efficiencies are given in Table |Vj All efficiencies in this 
article include constituent branching fractions. The yield uncertainties are statistical and 
uncorrelated systematic, respectively. The latter arise from modeling of multiple candidates 
in simulation and variations in the signal lineshape. Correlated systematic uncertainties are 
discussed separately in Section [V] 



TABLE V. ST yields and efficiencies including constituent branching fractions. Yield uncertainties 
are statistical and uncorrelated systematic, respectively, and efficiency uncertainties are statistical 
only. 



Mode 


Yield 


Efficiency (%) 


K--K+ 


75472 ± 300 ± 26 


63.74 ± 0.03 


K+TT- 


75655 ± 299 ± 26 


64.76 ± 0.03 


K+K- 


13813 ± 134 ± 5 


56.15 ±0.07 


7T + 7r~ 


6158 ±114 ±9 


72.08 ±0.11 


K° s ir°ir 


9209 ± 172 ± 16 


14.34 ± 0.04 


K° s tt 


23025 ± 174 ± 17 


31.53 ±0.04 


3251 ±81 ± 17 


10.81 ±0.05 


K s u; 


9292 ± 105 ± 7 


12.89 ±0.03 



B. Fully Reconstructed Hadronic Double Tags 

We reconstuct two categories of DT final states: 24 CP-allowed combinations of the 8 ST 
modes, and 136 modes with Kg7T + 7T~ . The second category includes 64 measurements where 
Yi and V, are paired with each of the eight ST modes [the sum of {7^~7r + , Y{\ and {K + tt~ , Y{\ 
(Cabibbo-favored), the sum of {K~n + ,Yj} and {K + n~ ,Yi} (Cabibbo-suppressed), and the 
sum of Yi and V paired with the six fully-reconstructed CP eigenstates], as well as 36 
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FIG. 2. ST M distributions and fits. Data are shown as points with error bars. The solid lines 
show the total fits, and the dashed lines show the background components. 



Cabibbo-favored combinations {YJ, Yj} + {%, Yj}, and 36 Cabibbo-suppressed combinations 
[Yi, Yj} + {Yi, Yj}. The 24 DT modes in the first category above were also used in Ref. [T2] , 
and we apply the same candidate selection criteria and yield determination methods as 
before, with the addition of the u mass sideband subtraction discussed above. Figure [3] shows 
some representative two-dimensional M distributions. Event counts in the signal regions (S) 
are corrected by background estimates from the sideband regions (A, B, C, D), with sideband 
scaling factors determined by fitting the ST distributions in simulated events in the same 
manner as Section HI A and integrating the fitted background function. Table VI gives the 
fully reconstructed DT yields and efficiencies for modes without KgTT + 7T~ . To account for 
systematic effects in the sideband definitions and in the extrapolation to the signal regions, 
we assign a 100% systematic uncertainty on the size of the sideband subtractions, which is 
much smaller than the statistical uncertainties in all cases. 



10 



For DT modes with K^-k+'k , we use the signal yields, efficiencies, and background esti- 
mates determined in Ref. [29] . 

3950912-003 
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FIG. 3. Two-dimensional M distributions with signal (S) and sideband (A, B, C, D) regions 
depicted, for {K~n + , K + ir~} and {K + K~, K^tt }. 




C. Double Tags with K£ 

For hadronic DT modes with a single K®, we employ the same partial reconstruction 
technique as for K^ir decays in Ref. [12], where the is identified by the four- vector 
recoiling against all other observed particles in the event. In the current analysis, we tag 
K^ii , K^T), K\uj, and K°7r°7r° decays with fully reconstructed ST candidates, as allowed 



by CP conservation and selected as described in Section HI A , and with an additional 
requirement of 1.86 GeV/c 2 < M < 1.87 GeV/c 2 . Each ST candidate is combined with a 
7T°, T), u candidate, or a pair of n° candidates. The signal process with K\ appears as a 
peak in the squared recoil mass, M^ iss , against this system. As in Ref. [12], we suppress 
the background by vetoing events with additional unassigned charged particles, but we veto 
addtional tt° candidates only for the K\r\ mode. In addition, for all K\ modes, we follow 
Ref. |29j by applying a veto on extra showers outside an energy-dependent cone around the 
predicted K\ direction. 

Figure shows examples of the resultant M^ iss distributions in data. We obtain yields 



from event counts in the signal and sideband regions as shown in Table VII, where the 
sideband is scaled by a factor determined from simulated events. We also subtract a small 
contribution due to continuum qq production, which is also estimated from simulated events. 



Table VIII lists the DT yields and efficiencies for K L modes without K s ix + ix . There 



are no uncorrelated systematic uncertainties for these modes. In the fit, we also include the 
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TABLE VI. Fully reconstructed DT yields and efficiencies including constituent branching frac- 
tions, for modes without K^tt^-k^ . Yield uncertainties are statistical and uncorrelated systematic, 
respectively, and efficiency uncertainties are statistical only. 



Mode 




Yield 


Efficiency (%) 


K-TT+ 


,K~7T + 


5.6 ±2.5 ±0.4 


41.5 ±2.8 


K-TT+ 


,K + TT- 


1731 ±42 ±11 


40.0 ±0.2 


K-TT+ 


,K + K- 


202 ± 14 ± 4 


35.2 ±0.5 


K-TT+ 


, 7T + 7T~ 


82.6 ±9.1 ±0.4 


44.5 ±0.9 


K--7T+ 


, K%TT°TT 


132 ± 12 ± 1 


8.6 ±0.2 


K--7T+ 


,K° s n° 


252 ± 16 ± 1 


19.4 ±0.3 


K--7T+ 


,K° s r) 


36.7 ± 6.2 ± 1.3 


6.9 ±0.3 


K--7T+ 


,K° s u; 


109 ± 11 ± 1 


8.5 ±0.2 


K+TT- 


,K+ir- 


4.0 ±2.0 ±0.0 


42.9 ±2.9 


K+ir- 


,K + K- 


191 ± 14 ± 1 


35.3 ±0.5 


K+ir- 


, 7T + 7T~ 


77.3 ±8.9 ±0.7 


45.6 ±0.9 


K+TT- 


,K° s TT°Tr 


121 ±11 ±2 


9.1 ±0.2 


K+TT- 


,K S 7T° 


242 ± 16 ± 


20.0 ±0.3 


K+TT- 




35.2 ±6.0 ±0.8 


6.9 ±0.3 


K+TT- 




89.4 ± 10.2 ± 1.3 


8.7 ±0.2 


K+K 




107 ± 11 ±2 


18.1 ± 0.5 


K+K 




24.6 ±5.0 ±0.4 


5.6 ±0.6 


K+K-,K° s uj 


47.6 ±7.2 ±0.0 


7.2 ±0.4 


7T + 7T~, 


K° s tt° 


37.0 ±6.1 ±0.0 


21.3 ±0.9 


7T + 7T~, 


KgT] 


6.0 ±2.5 ±0.0 


6.6 ±1.0 


7T + 7T~, 


K s u 


19.0 ±4.7 ±0.0 


9.4 ±0.7 


^7r°7r ,^7r° 


53.0 ±7.3 ±0.0 


4.1 ±0.2 




10.0 ±3.2 ±0.0 


1.4 ±0.2 


K° s 7r 7r ,K° s io 


18.0 ± 4.8 ± 1.0 


1.5 ±0.1 



TABLE VII. Signal and sideband regions in 

M miss (GeV 2 /c 4 ) for K° L modes. 



Mode 


Signal Region 


Sideband Region 


K° l tt° 


[ 0.10, 0.50 ] 


[ 0.80, 2.00 ] 


Kh 


[ 0.10, 0.45 ] 


[ 0.75, 1.75 ] 




[ 0.15, 0.40 ] 


[ 0.70, 1.25 ] 


K L 7T TT 


[ 0.10, 0.50 ] 


[ 0.80, 1.80 ] 



{Kg7r + 7c ,Kl7t } measurements from Ref. [29], which uses the same technique and criteria 
for Kin as those described above. 
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FIG. 4. Distributions of M^ iss for {if-7r+, i^vr } and {i^Tr , i^w}. The dashed lines show the 



signal and sideband regions. As indicated in Table VII the sideband for {K 7r + ,isT£7r } extends 
beyond the upper limit of the graph. 



D. Semileptonic Double Tags 



This Section describes our reconstruction of semileptonic D decays paired with fully 
reconstructed hadronic tags. In Section HI E below, we discuss an additional semileptonic 
mode with two undetected particles. 

In our previous analysis of Ref. [12], we reconstructed semielectronic final states inclu- 
sively, by identifying only the electron and not the accompanying neutrino or hadronic 
system. Also, we did not reconstruct semimuonic D° decays in Ref. [12]. In the current 
analysis, we replace inclusive reconstruction by exclusive reconstruction of K~e + u e and 
K + e~v e , which allows us, in general, to reduce systematic uncertainties because of lower 
background while keeping roughly the same statistical power. We also reconstruct K~ n + u^ 
and K + fi^u II without using the CLEO muon chambers because they are insensitive to muons 
in the momentum range of interest (p^ < 1 GeV/c). For both Kev e and Kfiv^, we begin 
with a fully reconstructed hadronic ST candidate, selected as described in Section [TTTA with 
additional requirements on M given in Sees. Ill D 1 and III D 2 below. Each ST candidate 
is then combined with a kaon and lepton candidate with opposite charges. To extract the 
signal yields, we calculate the quantity U = E m i SS — cp m i SS , where -E m i SS and p m i SS are the 
missing energy and the magnitude of the missing momentum, respectively, recoiling against 
the observed signal candidate particles in each event. Signal events peak at U — GeV 
because of the undetected neutrino. Representative U distributions for Kev e and K^iv^ are 
shown in Figs. [5] and [6j 
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TABLE VIII. DT yields and efficiencies including constituent branching fractions, for modes 
without KgTr + Tr~. Uncertainties are statistical only. All systematic uncertainties are correlated 

and are discussed in Sectio n PVT. 

Mode Yield Efficiency 



K-tt + ,KItt° 425 ±21 29.9 ± 0.2 

K + tt-,K° l tt° 381 ±20 31.1 ±0.2 

K%ir° , K° l tt° 235 ±15 14.6 ±0.1 

K° s 7],Kl^ 28.0 ±5.4 5.34 ±0.08 

K° s uj, K° l -k q 60.8 ±8.7 5.33 ± 0.07 

K-<K+,Klri 70.8 ±8.6 11.2 ± 0.1 

K + 7r~,Klrj 53.7 ± 7.6 11.5 ± 0.1 

K^TT°,K° L r] 21.7 ±4.8 5.55 ± 0.08 

K%r),K\r) 7.6 ± 2.8 2.01 ± 0.05 

K Q s u),Kl>r] 9.3 ±3.5 2.04 ±0.05 

K-ir + ,K° L u 143 ±13 12.1 ±0.1 

K+K-, K L u 155 ±14 12.2 ±0.1 

K^tt°,K l lo 80.7 ±9.8 5.70 ± 0.08 

K%ti,K%u 5.9 ±3.2 2.06 ±0.05 

K° s oj, K\uj 27.5 ±5.6 1.86 ±0.05 

K-ir+,Klir°iT 157 ± 13 13.0 ±0.1 

K + tt-,K1tt°tt 133 ±12 13.2 ± 0.1 

K + K~,KI-k -k q 57. 1 ±7.7 10.9 ±0.1 

vr + 7r-,K£7r 7r° 14.3 ± 4.9 14.5 ± 0.1 

K%Tr°ir°, Klir°ir° 36.6 ± 6.5 2.85 ± 0.06 



When the ST mode is K^ir , we form both the Cabibbo-favored (CF) combinations 
{K~tt + , K + £~i?g} and {K + tt~, K~£ + b>f} as well as the doubly Cabibbo-suppressed (DCS) 
combinations {K~ir + , K~£ + Vi} and {K + ir~, K + £~i)i}. As shown in Section [TT], the ratio of 
CF and DCS yields gives a direct measurement of the squared amplitude ratio r 2 . For the 
DCS combinations, we place an additional requirement on the polar angle of the kaon in the 
K~^7r^~ candidate of | cos 6 1 1 < 0.8. This requirement selects only kaons in the acceptance of the 
Ring Imaging Cerenkov counter (RICH), and it reduces the otherwise dominant background 
from misidentified Cabibbo-favored combinations to a negligible level. Figure [7] shows the 
summed U distributions for the two DCS Kev e modes and the two DCS Kfiu^ modes. 

For K£i/£ tagged with i^°7r + 7r~, we consider sums of CF combinations {%, K~l + v?} and 
{Yi, K + £~i'e}, as well as sums of Cabibbo-suppressed (CS) combinations {Yi, K~£ + ve} and 
{Yi, K + £~i'e}. We include the CF and CS {Yi,Kev e } yields measured in Ref. [29], and we 
perform a similar determination of CF and CS {Yi, K^u^} yields where we select Kg7i + n^ 
tags using the same criteria as in Ref. [29]. In analogy with tags, the ratios of CF 

and CS yields probes the amplitude ratios pf, integrated over each phase space bin. 
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FIG. 5. Distributions of U for {K~ir + , K + e~i? e } and {K%ir°, K~e + v e }. Data are shown as points 
with error bars. The solid lines show the total fits, and the dashed lines show the background 
components. 



1. Kev e 



For Kev e modes, we require the fully reconstructed ST to have 1.8530 GeV/c 2 < M < 
1.8780 GeV/c 2 for K° s ir° and K° s 7i 7i°, 1.86 GeV/c 2 < M < 1.87 GeV/c 2 for the DCS 
modes, and 1.8585 GeV/c 2 < M < 1.8775 GeV/c 2 for all other modes. We identify electron 
candidates with the same criteria as in Ref. [12], via a multivariate discriminant that com- 
bines information from the ratio of the energy deposited in the calorimeter to the measured 
track momentum (E/p), ionization energy loss in the tracking chambers (dE/dx), and the 
RICH. We fit the U distributions to signal and background probability distribution functions 
(PDFs) derived from simulated events and from sidebands in M and AE for the hadronic 
tag. For DCS modes, because of the low background, we simply count the number of events 
in the signal region \U\ < 0.352 GeV, and we estimate a relative background contribution 
of O(l0~ 3 ) from simulated events. 



Table IX gives the semileptonic DT yields and efficiencies for modes without Kg7r + 7T~ . 
The uncorrelated systematic uncertainties are determined from yield excursions under vari- 
ation of the signal and background shapes, histogram binning, and fit ranges. For DCS 
yields, we vary the size of the signal region by one-third, and we assign additional system- 
atic uncertainties to the kaon polar-angle requirement and electron identification. Yields for 
Kev e tagged with Kgii + TT~ (in 8 phase bins) are taken from Ref. [29] . 
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FIG. 6. Distributions of U for {K~ir + , K+pTv^} and {K%ir°, K'ii + u^}. Data are shown as points 
with error bars. The solid lines show the total fits, and the dashed lines show the background 
components. 



2. Kfx^ 

For all Kfiu^ modes, we require 1.86 GeV/c 2 < M < 1.87 GeV/c 2 for the fully re- 
constructed ST. We select muon candidates using the same criteria as for charged pion 
candidates, except the particle identification requirements (on dE/dx and RICH) are ap- 
plied to the muon mass hypothesis instead of the pion mass hypothesis. In reconstructing 
KfiUn, we reduce contamination from Kev e by requiring the muon momentum to be greater 
than 220 MeV/c, and we veto muon and kaon candiates that also satisfy the multivariate 



electron discriminant described in Section HID 1 Requiring p miss > 100 MeV/c suppresses 



the D — >■ Kir background, for which p m i SS peaks near zero. Finally, we reject events that 
contain an additional shower with energy greater than 100 MeV, in order to reduce the 
K~-k + -k° background contribution. After this requirement, K~tt + 7t remains the dominant 
background, but it is kinematically separated in U from signal K~[i + i>^ because of both the 
v^-it mass difference and the fi + -7i + mass difference. 

For all Kfiu^ modes, including the DCS modes, we determine signal yields by fitting 
the U distributions to signal and background PDFs derived from simulated events. For the 
non-DCS modes, the background PDFs are smoothed to reduce the effect of statistical fluc- 



tuations in the simulated histograms. Tables [X] and [XT] give the semileptonic DT yields and 
efficiencies for modes with and without Kg7r + n~, respectively. The uncorrelated systematic 
uncertainties are determined from yield excursions under variation of fit variations. For DCS 
yields, we also include systematic uncertainties for the kaon polar-angle requirement and the 
electron veto. 
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TABLE IX. DT yields and efficiencies including constituent branching fractions, for Kev e modes 
without Kgir + ir~ . Yield uncertainties are statistical and uncorrelated systematic, respectively, and 
efficiency uncertainties are statistical only. 

Mode Yield Efficiency (%) 

K-n + ,K + e-u e 1523 ± 40 ± 16 37.9 ± 0.2 

K+ir-,K + e-v e 5.0 ± 2.2 ± 0.9 30.6 ± 0.2 

K-K+,K+e~u e 156 ±13 ±4 33.0 ± 0.5 

iT-ir + ,K + e-ij e 70±9±5 42.3 ± 0.9 

K%TT°ir ,K + e-i> e 97 ±11 ±6 10.6 ±0.2 

K%^,K + eTv e 245 ±16 ± 15 20.1 ±0.3 

K° s r],K + e-u e 60 ± 8 ± 8 6.7 ±0.3 

K%u,K+e-v e 76 ±11 ±9 7.9 ± 0.2 

K-ir+,K-e+v e 9.0 ±3.0 ± 1.6 29.8 ± 0.2 

K+ir-,K-e + v e 1603 ± 42 ± 23 38.0 ± 0.2 

K-R + ,K-e + u e 175 ±14 ±8 33.6 ± 0.5 

7T-7T+, R-e + u e 64±8±1 42.2 ± 0.9 

K^TT°TT°,K~e + iy e 108 ± 12 ± 6 9.8 ±0.2 

K%Tr°,K-e+is e 244 ±16 ±5 20.2 ± 0.3 

K%n,K~e+v e 35±6±2 6.9 ± 0.3 

K° s uj, K~e + v e 73 ± 10 ± 8 7.5 ± 0.2 



TABLE X. DT yields and efficiencies including constituent branching fractions, for Kfiv^ modes 
without Kgir + ir~ . Yield uncertainties are statistical and uncorrelated systematic, respectively, and 
efficiency uncertainties are statistical only. 



Mode Yield Efficiency (%) 

K-ir + ,K+ii-v^ 1442 ± 40 ± 13 37.3 ± 0.2 

K+iT- ,K + fi~i? fl 7.0 ±2.7 ±1.1 34.8 ±0.2 

K-K + ,K + ^-v^ 121±12±0 32.9 ±0.2 

n-TV + ,K + n~i?^ 63.3 ±8.5 ± 1.0 42.7 ±0.2 

KgvrV, K+ii-Vp 85.2 ± 10.6 ± 4.7 8.6 ± 0.1 

K° s ^°,K + n~i>^ 216 ±16 ±6 18.3 ±0.1 

K^^K+ii-Dn 37.7 ±6.4 ±0.2 6.5 ±0.1 

K° s u,K + fi-9^ 91.9 ± 10.5 ± 1.6 7.1 ±0.1 

K-ir + ,K-fi+v^ 9.8 ±3.5 ± 1.6 33.8 ± 0.2 

K + TT-,K-fi+^ 1446 ±41 ± 13 38.0 ± 0.2 

K-K+jK-p+Vp 175±14±0 32.5 ± 0.2 

vr-vr + ,if^ + ^ 74.5 ±9.0 ±1.2 41.7 ±0.2 

K° s 7r°TT , K-fi+v^ 88.0 ± 10.5 ± 4.8 8.5 ± 0.1 

KgTr®, K~ fx^f^ 223 ±16 ±6 18.1 ±0.1 

K%r),K-n + Vn 33.0 ±6.2 ±0.2 6.5 ±0.1 

K%u,K-n+v^ 79.8 ± 10.3 ± 1.4 7.2 ±0.1 
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FIG. 7. Summed U distributions for {K ir + ,K e + f e } and {K + n ,K + e z/ e }, as well as for 
{iir~7r + ,if~At + ^} and {K + Tr~,K + n~u,j}. Data are shown as points with error bars. For Kev e , 
the vertical lines mark the signal region, the solid histogram shows the simulated signal distribu- 
tion normalized within the signal region to the number of observed events, and the small dashed 
histogram at large U values is the predicted background from simulation. For K/ii^, the solid 
histogram shows the total fit, and the dashed histogram shows the background component. 



E. Events with Keu e and K^tt 



To reconstruct {Kev e , K®ti }, where both the v e and K\ are undetected, we adopt the 
technique described in Refs. [37] and |38j for identifying events with two missing particles. 
Knowing the energy and momentum magnitude of the two D mesons in the inital state, and 
having reconstructed the K , e T , and 7r° in the signal process, the direction of each D meson 
in the e + e~ center of mass frame can be constrained to a cone around the flight direction of 
the K ± e T system or the ir°. In signal events, the D° and D° are collinear (within detector 
resolution), so when one of the cones is reflected through the origin, the two cones intersect. 
Background events typically have non-intersecting cones. 

Following Ref. [38J, we calculate the quantity 

x 2 D = 1 w— (cos 2 6 D Ke + cos 2 6 D T o + 2 cos 6 Ke cos 6 D Ke cos 6 D n o] , (16) 

sin 9 Ke ^o v ' ' 

where 6 Ke ^o is the angle between the K ± e T system and the 7r° candidate, and 6 D Ke and 6 D ^o 
are the opening half-angles of the D cones around the K ± e T system and the ir°, respectively. 
By construction, x 2 D is less than or equal to 1, and when the cones do not intersect, then 
x 2 D < 0. Signal events lie mostly in the range < x 2 D < 1, with a small tail extending to 
x 2 D < due to mismeasurement. 

In addition to the previously described criteria for K ± , e T , and 7r° candidates, we recover 
electron bremsstrahlung photons by augmenting the electron four-momentum by any showers 
located within 100 mrad of the track direction and unassigned to any other particle. Also, 
kaons that satisfy the electron identification requirements are rejected. 

We apply two additional requirements to suppress the background from radiative e + e~ 
Bhabha scattering events: the electron momentum is required to be less than 1 GeV/c, and 
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TABLE XI. DT yields and efficiencies including constituent branching fractions, for K/jlv^ modes 
with K^tt+it-. CF refers to the sum of Cabibbo-favored combinations {Yi,K and 
{Yi, K + fi~9^}. Similarly, CS refers to the sum of Cabibbo-suppressed combinations {Yi, K~ 
and {Yi, K + ^"9^}. Yield uncertainties are statistical and uncorrelated systematic, respectively, 
and efficiency uncertainties are statistical only. 



Mode Yield Efficiency (%) 



Y Q ,Knv^ 


(CF) 


162 ± 14 ±2 


18 


2 


±0 


2 


Y 1 ,K^iu fl 


(CF) 


75.7 ±9.3 ±0.8 


18 


4 


±0 


4 


Y 2 ,Kjiv l _ l 


(CF) 


132 ± 13 ± 1 


18 


6 


±0 


3 


Y^Kfii/f, 


(CF) 


36.3 ±6.4 ±0.4 


18 


9 


±0 


5 




{VjE ) 


f?7 7 4- 8 a 4- n 7 
1 . 1 ± o.o ± U. 1 


1 8 


y 


4- n 
± u 


q 




(CF) 


92.3 ±10.4 ±0.9 


17 


8 


±0 


3 


Y 6 ,Kfj,i/f, 


(CF) 


120 ± 12 ± 1 


17 


5 


±0 


3 


Y 7 ,Kfj,Vf, 


(CF) 


144 ± 13 ± 1 


17 


8 


±0 


3 


Yq,K^^ 


(CS) 


66.0 ±8.5 ±0.7 


17 


8 


±0 


4 


Y 1 ,Kfiu fl 


(CS) 


13.2 ±4.1 ±0.1 


17 


9 


±0 


6 


Y 2 ,K^,v^ 


(CS) 


33.0 ±6.2 ±0.3 


20 


9 


±0 


8 


Yz,Kiiv^ 


(CS) 


22.7 ±4.9 ±0.2 


20 





± 1 


1 


Y^KfiVfj, 


(CS) 


46.8 ±7.3 ±0.5 


18 


6 


±0 


6 




(CS) 


26.1 ±5.7 ±0.3 


18 


8 


±0 


7 


Yq,K^lv^ 


(CS) 


21.1 ±5.0 ±0.2 


17 





±0 


7 


Y 7 ,K\iv il 


(CS) 


58.1 ±8.2 ±0.6 


17 


7 


±0 


4 



the 7r candidate must not have daughter showers that lie within 100 mrad of the electron 
candidate. To reduce the background from {Kev e , i^vr } events, we reject events with extra 
tracks or TT°/r] candidates. 

Figure [8] shows the x 2 D distribution, which is fitted to signal and background shapes 
from simulated events. The fitted yield is 764 ± 36 ± 23, where the uncorrelated systematic 
uncertainty receives dominant contributions from the radiative Bhabha background deter- 
mination, variations in background and signal shapes, and uncertainties in the modeling of 
extra tracks and 7[°/r] candidates. The signal efficiency determined from simulated events 
is (34.58 ±0.04)%. 



F. Crossfeed and Peaking Backgrounds 



As in Ref . [12] , crossfeed among signal modes and peaking background contributions are 
subtracted by the fitter, using crossfeed probabilities and background efficiencies determined 
from simulated events, as well as branching fractions for peaking background processes 



These inputs are listed in Tables XII and XIII , and we account for correlated uncertainties 
among them. Backgrounds in DT modes are assumed to contribute at the same rate as for 
ST modes, with corrections for quantum correlation effects. When a ST background process 
does not exactly match the particle content of the corresponding signal process, then that 
background does not contribute to DT modes because of the additional constriaints imposed 
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FIG. 8. Distribution of x 2 D for {Kev e , K®ir }. Data are shown as points with error bars. The solid 
line shows the total fit, and the dashed line shows the sum of background components. 



by DT reconstruction. In general, peaking backgrounds contribute less than 1% to the 
measured yields, except in A"°7r° modes (1-2%), K\r\ modes (2-5%), \Kev e , K^tt } (1%), 
ST K%*°ir (3%), Cabibbo-favored {A^Tr+vr", K^tt"} (1-19%), and Cabibbo-suppressed 



{A>~ 



7T 



A°tT 



71 



} (2-56%). The large background fractions in the last two categories 



occur in modes with small numbers of observed events. 



IV. EXTERNAL MEASUREMENTS 



Unlike our previous analysis in Ref. [12], we do not include any external branching frac- 
tion or -Rws measurements in the fit. With our increased data sample, branching fraction 
measurements from other experiments do not have a significant impact on our precision, 
and our direct measurement of r 2 (from DCS {Kir, Kiu^}) obviates the need for input on 
r 2 from i?ws- We do, however, include the external measurements of mixing parameters 



shown in Table XIV, where x' 2 = 2Ru — y' 2 - Averages from Ref. [27] in this Table do not 



include previous CLEO-c results. Correlation coefficients that are not shown in Table XIV 



are taken to be zero. Also, since we neglect CP violation, we assume yep = y- In the fit, 
we use the signed measurement of x along with y' to resolve the sign ambiguity in sin 5 and 
the Si, even though we quote a fitted value only for x 2 . 
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TABLE XII. Crossfeed probabilities among signal modes. Ranges are given for groups of modes. 



Crossfeed — > Signal 



Efficiency (%) 



TT 



K-TT+ 



K^ 



o 



,K° sV 
,K° sV 



K~7T+ 

K+tt~ i 

K~TT+ 
K+TT~ 
K-TT+, 
K+TT- 
K° S 7T , K° l tt° 
K° s v, K° l tt° 

-0 

K-ir 

K+TT 

K--K 

K+TT 
K-TT 
K+TT 

K+K 



TT ' TT 
K° S TT°TT 
K-TT+ 
K+TT- 

Kevp 



K-TT+ 
K+TT- 
K-TT+ 
K+TT~ 

K^ 
K~ 

K+TT- 



TT 
'TT' 
"7T H 



K° L TT° 



K° r TT° 



K° s uj 
K° s u 

K° S TT°TT 
K° S TT°TT 
K° S TT° 
K° S TT° 
K° S TT° 
K° S TT° 
K° S TT° 
K° L TT° 
K° L TT° 

~e v e 
~e+u e 
'e~u e 
~e+v e 



K 
K 
K 
K~ 
K+ 
K- 
K^ 
K- 



K^rj, K/iV/j, 
K^uj, K/iv^ 



K~TT+ 
K+TT- 
K S TT°TT° 
K'TT 
K+TT 
K-TT 

K+ 



TT 



TT 



K° L V 

Kin 

K° L rj 
K° L rj 
K° s tt°, Kl V 
K° sV , K° lV 

K° L V 



K 



K+TT- 



K° L TT° 
K° L TT° 



K-TT+ 
K+TT- 
K-TT+ 
K+TT- 
K-TT+ 
K+TT- 

K+K- 

TT+TT- , 



K° S TT°TT , 
K-TT 



K+TT- 



Keu e , 

K-TT+ 
K-TT+ 
K+TT- 

K+- 



K° L u 
K° L u 

K° L TT°TT 
K° L TT°TT 
K° L TT°TT 
K° L TT°TT 
K° L TT°TT 
K° L TT°TT 
K° L TT ^ 
K° L TT°TT 
K° L TT°TT 
K°TT° 



TT 



TT 



K-TT+ 
K+TT- 
K+TT- 



K' 
K~ 
K' 
K' 

K- 



e 
e 
e 

e~v t 
\i+v. 



K-jj+v, 



K- 



K° S TT- 



Kfj,Ufj, 
Kjiv^ 



0.088 ± 0.002 
0.089 ± 0.002 
0.080 ± 0.003 
0.45 ±0.02 
0.43 ±0.02 
0.13 ±0.01 
0.12 ±0.01 
0.14 ±0.01 
0.13 ±0.01 
0.06 ±0.01 
0.04 ±0.01 
0.03 ±0.01 
0.13 ±0.01 
0.12 ±0.01 
0.47 ±0.02 
0.46 ±0.02 
0.61 ±0.03 
0.62 ±0.03 
0.49 ±0.02 
0.68 ±0.03 
0.16 ±0.01 
0.11 ±0.01 
0.12 ±0.01 
2.31 ±0.03 
0.0048 ± 0.0022 
0.0057 ± 0.0025 
0.0038 ± 0.0020 
0.0019 ±0.0014 
0.0029 ± 0.0017 
0.0067 ± 0.0026 
0.0019 ± 0.0014 
0.0010 ±0.0010 
(0-3) x 10~ 3 
(0-5) x 10~ 2 



V. SYSTEMATIC UNCERTAINTIES 



We include systematic uncertainties directly in the fit. Uncorrelated uncertainties for 



each yield are discussed above in Section III Correlated uncertainties are given in Tables XV 



and |XVI| When the fit is performed without external measurements, these correlated uncer- 
tainties cancel in all the fit parameters except M and the branching fractions. Uncertainties 



on the peaking background branching fractions in Table XIII and the external measurements 
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TABLE XIII. Peaking background branching fractions and efficiencies including constituent branch- 
ing fractions. Ranges are given for groups of modes. Fully reconstructed hadronic tag modes are 
denoted by X. Backgrounds marked by an asterisk (*) occur only in STs, not DTs. 



Background — > Signal 



Sbkg (%) M Efficiency (%) 



D+ 



7T 

*° (*) 
*° (*) 

A° 

Generic D°D° (*) 
Generic D + D~ (*) 
r]-K°, X 



7]TT°, 



Keu, 
K°tt 



7T 



7T 



X 
X 
X 

K° L TT°7T 
7T - , X 



K° S 7T°7T 
K S 7T 7T 
K S 7T°7T 
K° S 7T° 
K° S 7T° 
K° S 7T° 
K S 7T° 
K° L TT° 
K° L 7T° 

K° L V, 

K° S TT+ 



, x 

X 
X 



X 



2.94 ±0.16 
13.9 ±0.5 
6.90 ±0.32 
1.447 ±0.046 
0.373 ±0.022 
100 
100 

0.064 ±0.011 
0.080 ±0.008 
0.064 ±0.011 
0.167 ±0.019 
From fitter 
2.94 ±0.16 
0.744 ±0.021 



0.0076 ± 0.005 
0.0027 ±0.0001 
0.0594 ±0.0010 
0.078 ±0.004 
0.011 ±0.004 
0.0006 ± 0.0001 
0.0003 ± 0.0002 
0.3-1.4 
0.8-4.3 
0.1-0.3 
0.2-0.8 
0.45 ± 0.04 
(0-6) x 10" 3 
(1-25) x 10~ 2 



TABLE XIV. External measurements of yep, x, y, r 2 y' and x' 2 . Correlation coefficients p are for 
the measurements in the same column. 



Parameter 




Value(s) (%) 




VCP 




1.064 ±0.209 [27] 




X 




0.419 ±0.211 [27] 




y 




0.456 ±0.186 [27] 






0.364 ±0.017 [8] 


0.303 ±0.016 ±0.010 [9] 


0.304 ± 0.055 [TO] 


y' 




0.97 ±0.44 ±0.31 [§] 


0.85 ± 0.76 HU] 


x' 2 


o.oislHia IB! 


-0.022 ±0.030 ±0.021 [9] 


-0.012 ±0.035 [E] 


p(r 2 ,y') 


-83.4 gO] 


-87 00] 


-97.1 gO] 


p(r 2 ,x' 2 ) 


±65.5 00] 


±77 g0] 


±92.3 gO] 


p(y',x' 2 ) 


-90.9 00] 


-94 gO] 


-98.4 gO] 



in Table XIV are also included as systematic uncertainties. 

In Table XV , we list the correlated systematic uncertainties on reconstruction and parti- 
cle identification efficiencies for individual final state particles. Those for tt , K , it , Kg, 
and e reconstruction and particle identification (PID) are determined using the techniques 
described in Refs. [U]II2]- ^ or we inciu de additional contributions for the flight signifi- 
cance and invariant mass requirements. For rj, we use the same uncertainties as in Ref. 



The K\ uncertainties are determined by varying the selection requirements. The above 
studies also provide efficiency corrections, which are included in the signal and background 
efficiencies listed in this paper. 

Table [XVI] shows mode-dependent uncertainties, including those for initial state radiation 
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TABLE XV. Correlated, fractional efficiency systematic uncertainties and the schemes for applying 
them in the fit. 



Source Uncertainty (%) Scheme 

Track finding 0.3 per track 

hadronic interactions 0.5 per 

K° s finding 0.9 per K Q S 

7T° finding 2.0 per tt° 

r) finding 4.0 per rj 

dE/dx and RICH 0.1 per vr ± PID cut 

dE/dx and RICH 0.1 per PID cut 

Electron identification 0.4 per 

shower veto 0.4 per K\ 

background subtraction 0.7 per K\ 

K° L track veto 0.3 per K\ 

signal shape 1.4 per K\ 



(ISR), final state radiation (FSR), and selection requirements in two-track modes that veto 
radiative Bhabhas and cosmic muons. Most of the uncertainties in Table IXVll are determined 
with the same techniques as in Ref. p2], updated to the current data and simulated samples. 
We assess uncertainties for new modes by relaxing selection requirements and noting the 
resultant change in efficiency-corrected yield. 



VI. FIT RESULTS 



We combine the 261 yield measurements discussed above, along with estimates of effi- 
ciencies and background contributions, in a \ 2 fit that determines 51 free parameters: M\ y; 
r 2 ; cos 5; sin 5; x 2 ; p 2 , Ci, and Sj for each of the 8 phase bins in Kg7c + 7C~; and 21 branching 
fractions. We tested the analysis technique using a simulated sample of quantum-correlated 
D°D° decays with an effective integrated luminosity of 10 times our data sample. In Ta- 



bles [XVTT and XVIII, we show the results of two fits: one with no external inputs (Standard 
Fit) and one including the measurements in Table XIV (Extended Fit). The fits are per- 



formed using both statistical and systematic uncertainties on the input measurements. We 
evaluate statistical uncertainties on the fit parameters by performing a second set of fits using 
only statistical uncertainties on the inputs. Then, we compute the systematic uncertainty 
on a given fit parameter by taking the quadrature difference between its total uncertainty 
and its statistical uncertainty. When a fit parameter is directly constrained by an exter- 
nal measurement, we quote only one uncertainty. Asymmetric uncertainties are determined 



from the likelihood scans discussed below. Table XIX gives the correlation coefficients for 



y, r 2 , cos 5, sin 5, and x 2 . The full correlation matrices can be found on EPAPS |43j . 

These results include the first direct measurement of sin 5, as well as first branching frac- 



and K° L 7i 7i° 



. We also present the first CLEO-c 
In the Standard Fit, the statistical 



tion measurements for D° — > K^rj, K^co, 
branching fraction measurement for D° — > K~ jj^ 
uncertainties on y and r cos 5 are approximately three times smaller than in our previous 
analysis [12]. Because of the strong correlation between r 2 and cos 5, and because r 2 is 
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TABLE XVI. Correlated, mode-dependent fractional systematic uncertainties in percent for STs. 
An asterisk (*) marks those uncertainties that are correlated among modes. The schemes by which 
these uncertainties are applied to DTs are also given, along with the fractional DT uncertainty, 
AdTj on mode {A, B} for ST uncertainties of a on mode A and /3 on mode B. 





AE 


ISR* 


FSR* 


Lepton Veto* 


Other 




0.5 


0.5 


0.9 


0.5 




K+K~ 


0.9 


0.5 


0.5 


0.4 


0.5 K± cos 9 cut 


7T + 7T~ 


1.9 


0.5 


1.4 


3.2 






2.6 


0.5 






1.5 Kg daughter PID 
0.7 resonant substructure 




0.9 


0.5 








K° s v 


5.5 


0.5 






0.3 jj mass cut 
0.7 B(ry -»• 77) [39] 




1.2 


0.5 


0.6 




0.1 w mass cut /SB subt. 
0.8 B{uj -4 7T + 7T"7r ) [39] 


K£vr°(7r ) 




0.5 








k° lV 




0.5 






1.6 extra 7r° veto 






0.5 


0.6 




0.1 w mass cut /SB subt. 
0.8 B(u ->• tt+tt-tt ) [39] 


K° S TT + TT- 


0.9 


0.5 


1.4 






Keu e 




0.5 


0.3 




2.0 spectrum extrapolation 






0.5 


0.3 




2.0 spectrum extrapolation 
0.4 extra shower veto 


Scheme 


per D 


per yield 


per D 


per ST 


per D 


Adt 


Va 2 + /3 2 


(a + /3)/2 


a + /3 





y/a 2 + P 2 



determined by different inputs in the two analyses, a direct comparison of the cos 5 uncer- 
tainties is not instructive. Accounting for correlated uncertainties, the results for cos 5 in 
the Standard and Extended Fits differ by 2.5a. This difference is caused primarily by the 
upward fluctuations in y and r 2 in the Standard Fit, which are both negatively correlated 
with cost5. 

The above factor of three improvement in y and r cos 5 can be attributed in equal parts to 
the increased size of the dataset and to the additional final states in the current analysis. In 
particular, the new K\ modes (including {Kev e , Kin }) reduce the statistical uncertainties 
on y and rcos5 in the Standard Fit by roughly 10%, and the Kfiu^ modes reduce them by 
20-30%. The K^ir + n~ measurements from Ref. [29] have a similar effect on y and rcos<5 as 
the new K® modes, but they also provide all the information on sin S in the Standard Fit. 

The Extended Fit demonstrates that when our cos 5 and sin 5 measurements are used 
to combine y and y' from other experiments, the overall uncertainty on y is reduced by 
approximately 10% compared to the global average of all measurements in Table XIV| (except 
yep from Ref. [5]) found in Ref. [10]: y = 0.79 ±0.13. Note that this global average includes 
the results from Ref. [12] . 

Figure [9] shows the one-dimensional posterior PDFs for cos S, sin 5, 5, and y in the Stan- 
dard Fit, including statistical and systematic uncertainties. These curves are obtained by 
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TABLE XVII. Results from the Standard Fit and the Extended Fit for all parameters except 
branching fractions. Uncertainties are statistical and systematic, respectively. In the Extended 
Fit, we quote only one uncertainty for y, r 2 , and x 2 , which are directly constrained by an external 
measurement. 



F*a T*a in ptpr 

J. CliA CAj 111 XJ v_> 1 




Extended Fit 

J i-/V Av K^j ± 1 v_> VI 1 lu 


Af (\0 6 ) 

J V I XVJ 1 


m . v > u — _i_ vj . vjo vj _i_ vj.vj^x-vj 


q 1 1 4 _i_ 050 + 040 

O . J. 1 1 _l_ VJ . VJtJ VJ _l_ VJ.VJ'XVJ 


v (%) 


A o _i_ o n _i_ i n 

J: . ■ I . VJ i X . yj 


fi'iK + 0114 

\ ) .yj '.) \> _i_ vj . x x 'x 


r 2 (%) 


n coo + (] 1f)7 + f) 045 

VJ.vJtJtJ _l_ U.1UI _l_ v7 ■ \J~tiJ 


QQQ + 008 

VJ . 'J (J 'J _l_ VJ . VJVJO 


OfVl n 
uUu u 


u - ol -0.18-0.05 


1 1 c+0.19+0.00 
± -- Ld -0.17-0.08 


cm A 

QllX u 


_n ni + o 41 + n D4 

VJ • VJ X I VJ . t: X I VJ . VJ t: 


u - tJU -0.31-0.20 


r 2 (%) 


n dfi + n 2^ + n 1 1 

VJ . VJ VJ I VJ . — • O I VJ . X X 


0022 + 0023 

VJ . \J\J £1 £1 1 VJ . VJ VJ ^1 


PO 


337 ± (130 ± D 006 

V> • i J ' J 1 _1_ v7 . VJ O VJ _1_ \J •\J\J\J 


352 ± 032 ± 005 

vj.otj^ _i_ v . y ) >._> ±-i _i_ vj.vjvjtj 


ol 
Pi 


n ovn _i_ n 044 + D 005 

\J.Li I VJ _l_ VJ . \J L ± L ± _1_ W.VJWtJ 


o«o + 047 + 000 

VJ.iiOVJ _l_ VJ.yJ'-k I _1_ VJ.VJVJVJ 


ol 

P2 


VJ.Zj*JtJ _l_ VJ . VJZjO _l_ V / . V 7 V 7 'J 


952 + 028 + 004 

VJ • ZjfJ Zj I VJ . VJi O I VJ . VJ VJt: 


P3 


n qoo + n 066 + n 005 


41 fi + 069 + 000 

V J . t: X VJ _l_ VJ . VJ VJ sJ _l_ V J . V> VJ V J 


O 2 , 
Pi 


592 ± 067 4- 010 

\J . tj iJ .-J _1_ VJ . VJ VJ 1 _1_ VJ . VJ X VJ 


623 ± 071 ± 000 

vj . vjzj *j _i_ vj . vj 1 x _i_ \j •xjyjyj 


2 r 

P5 


U.U^O _l_ VJ.VJ'X-'X- _l_ L/.VJVJL/ 


ooq + 040 + 008 

V / . '..) >"J _l_ VJ . VJt: VJ _l_ VJ. VJVJO 


P6 


nufi + n 023 + o ooo 


1 45 + 023 + 000 

\J . 1 t: '-.J _l_ V7 . \J Z. k) _l_ VJ.VJVJVJ 


in 


445 + 039 + 002 

\J . 'J _1_ VJ.VJvJtJ _l_ VJ . VJ VJ.£j 


4-10 + 039 + 003 

V . ttO >'J _l_ VJ.VJOiJ _l_ V . VJVJ tJ 


' 


—0 76 ± 06 ± 01 

VJ . 1 _1_ VJ.VJVJ _1_ \ I .\f 1 


— 73 ± 06 ± 01 

\> • 1 ' ) _1_ VJ . VJVJ —1— VJ . VJ X 


'1 


_0 75 ± 1 1 + 00 

VJ . 1 <J _l_ U.ll _l_ VJ . VJVJ 


_0 79 + 1 1 + 02 

VJ . 1 _ I VJ . X X I \J •\J j—i 


On 


00(1 i ninxom 

VJ . VJ VJ _l_ VJ . 1 VJ _1_ \J . V 7 1 


03 + 1 + 02 

VJ.VJtJ _l_ VJ . X VJ _l_ VJ.VJX 




VJ.^(J _l_ VJ.1(J _l_ VJ . VJ X 


47 + 1 4 + 01 

VJ . t: 1 I VJ . X t: _1_ VJ . VJ X 




05 + o 07 + 01 

VJ.tJtJ _l_ VJ . VJ i _l_ VJ.VJX 


05 + 07 + 00 

VJ.tJtJ _l_ VJ . VJ 1 _l_ VJ . VJ VJ 


< 5 


70 + 0Q + 01 

Vl . 1 iJ _l_ VJ . VJ tj _l_ VJ.VJX 


»i + 09 + 00 

VJ.OX _l_ VJ.VJtJ _l_ VJ . VJ VJ 


Cfi 
Hb 


-0.20 ± 0.13 ± 0.02 


-0.16 ± 0.13 ± 0.01 


C7 


-0.41 ±0.07 ±0.01 


-0.39 ±0.07 ±0.01 


so 


0.55 ±0.16 ±0.00 


0.61 ±0.15 ±0.02 


si 


0.53 ±0.28 ±0.00 


0.56 ±0.27 ±0.03 


•52 


0.93 ±0.15 ±0.00 


0.91 ±0.15 ±0.02 


S3 


0.47 ±0.30 ±0.00 


0.52 ±0.29 ±0.01 


S 4 


0.55 ±0.24 ±0.00 


0.60 ±0.23 ±0.02 


S5 


-0.71 ±0.24 ±0.00 


-0.69 ±0.24 ±0.00 


S6 


-0.42 ±0.27 ±0.06 


-0.17 ±0.29 ±0.03 


«7 


-0.30 ±0.18 ±0.04 


-0.21 ±0.19 ±0.03 


Xfit/ndof 


193.2/210 


214.7/222 



re-minimizing the x 2 & t each point and computing £ = e~( x2 ~ x ™^^ 2 . Also shown are the 
two-dimensional contours for combinations of y, cos 5, and sin 5. Because of the sin 8 sign 
ambiguity, the PDFs for both sin S and 8 in the Standard Fit are symmetric around zero. 
Figure 10 shows the same distributions for the Extended Fit. Here, the sign ambiguity is 
resolved by the external measurements. All the PDFs, except those for 8, are well described 
by Gaussians or bifurcated Gaussians. In particular, the non-Gaussian profile for cos 5 in 
the Standard Fit from Ref. [12] has been eliminated by our direct measurement of r 2 . 
Although the central value for cos 5 in the Extended Fit is unphysical, we find r 
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TABLE XVIII. Branching fraction results from the Standard Fit and the Extended Fit. Uncer- 
tainties are statistical and systematic, respectively. 



Parameter 


Standard Fit 


Extended Fit 


B(K-tt + ) (%) 


3.77 


±0.06 


± 


0.05 


3.76 


± 


0.06 


± 


0.05 


B(K-K+) (1(T 3 ) 


3.99 


±0.07 


± 


0.08 


3.98 


± 


0.07 


± 


0.08 


B(ir-ir + ) (1(T 3 ) 


1.36 


±0.03 


± 


0.04 


1.37 


± 


0.03 


± 


0.04 


B(K° S TT°TT ) (%) 


0.99 


±0.02 


± 


0.06 


0.99 


± 


0.02 


± 


0.06 


B{Kl^) (%) 


0.94 


±0.03 


± 


0.03 


0.96 


± 


0.03 


± 


0.03 


BiKU (10- 3 ) 


3.36 


±0.30 


± 


0.17 


3.40 


± 


0.31 


± 


0.17 


B(K° L u) (%) 
B(K° S *°) (%) 


0.90 


±0.05 


± 


0.03 


0.91 


± 


0.05 


± 


0.03 


1.17 


±0.02 


± 


0.03 


1.16 


± 


0.02 


± 


0.03 


B(K° s r,) (10- 3 ) 


4.95 


±0.14 


± 


0.36 


4.90 


± 


0.14 


± 


0.36 


B(K®ui) (%) 


1.15 


±0.02 


± 


0.04 


1.14 


± 


0.02 


± 


0.04 


B(Klir°TT ) (%) 


0.95 


±0.06 


± 


0.05 


0.94 


± 


0.06 


± 


0.05 


B{K~e + v e ) (%) 


3.54 


±0.05 


± 


0.08 


3.52 


± 


0.05 


± 


0.08 


B(K-fi+u^) (%) 


3.38 


±0.05 


± 


0.08 


3.36 


± 


0.05 


± 


0.08 


H(F ) (10- 3 ) 


4.38 


±0.18 


± 


0.12 


4.33 


± 


0.17 


± 


0.11 


B(n) (10- 3 ) 


1.65 


±0.10 


± 


0.04 


1.63 


± 


0.10 


± 


0.04 


s(y 2 ) (10- 3 ) 


3.43 


±0.16 


± 


0.10 


3.33 


± 


0.14 


± 


0.08 


fi(y 3 ) (10- 3 ) 


0.99 


±0.08 


± 


0.03 


0.97 


± 


0.08 


± 


0.02 


B(Y A ) (10- 3 ) 


1.70 


±0.11 


± 


0.05 


1.62 


± 


0.10 


± 


0.04 


s(y 5 ) (10- 3 ) 


2.11 


±0.13 


± 


0.07 


2.13 


± 


0.12 


± 


0.05 


fi(y 6 ) (10- 3 ) 


3.15 


±0.15 


± 


0.08 


3.14 


± 


0.14 


± 


0.08 


B(Y 7 ) (10- 3 ) 


3.68 


±0.16 


± 


0.09 


3.71 


± 


0.16 


± 


0.09 



(cos 2 5 + sin 2 5) 1 / 2 = 1.28 ± 0.27 to be consistent with physical boundary. Similarly, in the 
Standard Fit, r = 0.81 ± 0.21. The PDFs for 5 in Figs. [9] and 10 are obtained by probing 
cos S and sin 5 under the constraint r = 1, which reduces the height of the PDF relative to 
the other parameters. The implied values for S from these PDFs are \S\ = (lOt 2 ^ 3 ) f° r the 
Standard Fit and S = (18^7)° for the Extended Fit. Also, applying the above constraint 
in the Standard Fit improves the uncertainties on y and x 2 by 15% and 8%, respectively, 
resulting in y = (3.3 ± 1.7 ± 0.8)% and x 2 = (0.14 ± 0.21 ± 0.09)%; the changes in all other 
parameters are negligible. Performing the Extended Fit with r = 1 produces negligible 
shifts in all the fit parameters. 

Our results for q, Sj, and branching fractions do not supersede other CLEO-c measure- 
ments. For Cj and Sj, our fitted values are consistent with those in Ref. after accounting 
for differences between the two analyses. 



VII. SUMMARY 

We present an updated analysis of quantum correlations in D°D° decays at the ^(3770) 
using the full CLEO-c dataset, resulting in a new value of cos 5 = 0.8lto^l-ox)l an d a first 
measurement of sin 5 = —0.01 ± 0.41 ± 0.04, which, when combined, imply a strong phase 
of \5\ = (lO^jljto 3 ) . By including external inputs on mixing parameters in the fit, we find 
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TABLE XIX. Correlation coefficients (%) for the fits in Table XVII using both statistical and 
systematic uncertainties. 



r 

cos 6 
sin 5 

x 2 

r 2 

cos 5 
sin 5 

x 2 



y 



cos 5 sin 5 



Standard Fit 


-53 -42 

-3 +1 +4 
-73 +39 +2 

Extended Fit 

+3 
-27 -16 
+62 +21 +36 

+9 +25 +5 -18 



alternative values of cos 6 = 1.15±g;i?±8;8g, sin 5 = 0.56±g;|f±g;lJ, and 5 = (18±#)°. The 
effect of these measurements on the world averages of mixing parameters is to improve the 
uncertainty on y by approximately 10%. 
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